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Abstract 
One of the most important theories in linear programming is the dualistic theory and its basic idea is that for every 
linear model has dual linear model, so that solving the original linear model gives a solution to the dual model. 
Therefore, when we solving the linear programming model, we actually obtain solutions for two linear models. In 
this research, we present a study of the models. The neutrosophic dual and the binary simplex algorithm, which 
works to find the optimal solution for the original and dual models at the same time. The importance of this algorithm 
is evident in that it is relied upon in several operations research topics, such as integer programming algorithms, 
some nonlinear programming algorithms, and sensitivity analysis in linear programming... 
Keywords: Neutrosophic Science; Neutrosophic Linear Models; Dual Models; Neutrosophic dual Linear Models; 
binary Simplex Algorithm. 
1. Introduction: 
Most companies and institutions rely on studies provided by experts and researchers using operations research 
methods in order to ensure a safe work environment away from danger and danger and to give decision makers a 
margin of freedom, after the studies and research presented using the concepts of neutrosophic science in most fields 
of science have proven their ability to provide more accurate results. From the results that we were obtaining using 
classical data-driven studies, we reformulated many operations research topics using neutrosophic concepts [1-20], 
and to complement what we have done, we present in this research a study of dual linear models and the binary 
simplex algorithm that gives us a solution for the original and dual models at the same time so that we can provide 
a clear study that helps decision makers in companies and institutions develop plans and programs through which 
the greatest profit and lowest cost are achieved. 
Discussion: 
In our practical life, we encounter many problems that are formulated in the form of linear mathematical models 
consisting of an objective function and a set of constraints in the form of equations or inequalities. The linear 
model is written in many formulas that are distinguished by the type of the objective function and the form of the 
constraints. Many references have been provided that studies of linear and dual linear models are based on classical 
data, and the appropriate algorithms for solving each of them and how to find the dual model for any model [21- 
25]. Since classical data results in optimal solutions, which are classical values, restricted values appropriate to the 
conditions in which the data were collected, and for any change in these conditions, these solutions are 
inappropriate and may cause the facility a large unexpected loss, so we found it appropriate to reformulate these 
studies using concepts neutrosophic science is done by taking data with neutrosophic values suitable for all 
conditions. In previous research, we presented neutrosophic linear models, the direct neutrosophic simplex 
algorithm, and the graphical method for finding the optimal solution for linear models,[18,19,20] . As a 
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continuation of what we have done previously, we present in this research a study whose purpose is to reformulate 
the binary simplex algorithm to find the optimal solution for the original model and the dual model at the same 
time. 


We present this study in two stages: 
The first stage: 
Finding the neutrosophic dual models using the binary table. 


When we want to find the model associated with any linear model, we start by placing this model in the 
symmetrical form, so we mention the symmetrical form. 


The symmetrical form of the neutrosophic linear model: [19] 


We say of a linear model that it is in the symmetrical form if all variables are constrained to be non-negative and 
if all constraints are given in the form of inequalities (and the inequalities of the maximization model constraints 
must be in the form (< less than or equal to) while the inequalities of the minimization model constraints must be 
in the form (2 is greater than or equal to), then the linear model is written in the neutrosophic symmetric form in 
one of two cases: 


First case: 


n 
NZ = YG + &)x; — Max 
i=1 


Constraints: 


Najyx; <b; 46; ; i=1,2,..,m 


n 


jJ=1 


Second case: 


j=1 
Constraints: 


In both cases, we have xj 20, which are the decision variables, unknown values that we obtain after solving the 


linear model. 


Ne; = c; + €; and Nb; = b; + 6; and Na;; = a;; + wy; where G=1,2,....n, i=1,2,...,m) are the data of the problem 


under study, and they are neutrosophic values, indefinite values that enjoy a margin of freedom and are taken 
according to The nature of the situation represented by the linear model 


First method: 
Constructing neutrosophic dual linear models using tables: 


If the original model is given in detailed form, we put it in the symmetrical form as stated in the previous paragraph, 
then we draw a binary table for the original and dual models according to the following steps: 
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1- The coefficients of the objective function in the original model are the constants column in the dual model, 
and the constants column in the original model are the coefficients of the objective function in the dual 
model. 

2- We invert the signs of the inequalities of the constraints (if they were in the original model of type < they 
become in the dual model of type =) 

3- We change the objective from maximizing in the original model to minimizing in the dual model 

4- We place each constraint (row) in the original model corresponding to a column in the dual model, 
meaning there is one variable for each constraint in the original model 

5- The variables in the original model and the dual model satisfy the non-negativity constraints. 

We explain the above through the following two cases: 


The first case: The original model is symmetrical and of the maximization type 
Find 
NZ = Neyx, + NegxX2 + +++ + NeyX, — Max 


Constants: 
Nay1X1 + Nay2X_ ++: +NQynX%y, < Nb, 


Naz4X4 + Nay2X2 tet NQznXn < Nb, 


NQyn1X%1 + NAm2X2 + + NAmn%y < Nbm 


X41) XQ, 00) Xn 20 


The binary table for the original model and the accompanying model is as follows: 
Table No. (1) Objective follower of the maximization type 


objective 


Nay4X4 + Nay2X2 ++. 
Naz4X 1 + Naz2X2 ++ 


m Nm X1 + NAm2X2 + +° + NAmnXn 


IV IAT IA} IA FIA 


Non-negativity 
constraints 


Kay Maire Ny 


Objective 


Nay1¥1 + NQgiy2 +++ + NAmiVm 


NQy2y1 + Nagoy2 +++ + NAm2Vm 


Nainyi + N@any2 ++ Nan 
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Non-negativity > 
constraints 


The second case: The original model is symmetrical and of the reduction type: 
Find 


NL = Neyx1, + NegX2 +++ + NCyXn — Min 
Constants: 
Nay4X 4 as Nay2X2 + ve + NQinXn > Nb, 


Nap4X4 + NaQp2X2 tet NQznXn => Nb, 


NQm1%1 + NAm2X2 +++ NAmnXn = Nbm 


X4,X2,+1,Xy, 20 


The binary table for the original model and the accompanying model is as follows: 
Table No. (2) objective follower in the original model of the reduce type 


objective 


a Nari +Nayt2 t+" tNGintn | > | 


Nayz4X1 + Naz2X2 +++: + Nag Xp 


IV IV} IV 


NQmX1 + Negetece + NAmnXn 


ee 


V 


Non-negativity 
constraints 


objective 


2 NQy2V1 + Nagzy2 +++ + NOm2Vn 


IA 


Nainyi + Naany2 ++ NdmnYm 


Non-negativity 
constraints 


The second stage: formulation of the binary neutrosophic algorithm. 


The binary simplex algorithm (for the original and dual models) neutrosophic. Through this algorithm, we can find 
the two ideal solutions for both the original and dual models at the same time. Before starting the binary simplex 
algorithm, we must mention the modified simplex algorithm that we will use within the steps of the binary 
algorithm. 
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Modified simplex algorithm: 


In the modified Simplex algorithm, after converting the regular linear model to the basic form, we place the 
coefficients in a short table whose first column includes the basic variables and whose top row includes the non- 
basic variables only. We define the pivot column, which is the column corresponding to the largest positive value 
in the objective function row if the objective function is a maximization function (but if the objective function is 
a minimization function, itis the column corresponding to the most negative values). Let this column be the column 
of the variable x,. We define the pivot row. The pivot row is determined. Through the following indicator: 


Nb;]__ Ny 
Na, ~ Nats 


No = min| > 0; Na;; > 0,Nb; > 0 

Let this line be the line of the base variable y;, then the anchor element is the element resulting from the intersection 
of the anchor column and the anchor line, i.e., Na,,. Then we calculate the new elements corresponding to the 
anchor line and the anchor column as follows: 


1- We put opposite the pivot element Na,, the reciprocal of — 
ts 


2- We calculate the elements of the row corresponding to the pivot row (except the pivot element) by 
dividing the elements of the pivot row by the pivot element Na;, 

3- We calculate all the elements of the column opposite the pivot (except the pivot element) by dividing the 
elements of the pivot column by the pivot element Na,, and then multiplying them by (-1) 

4- We calculate the other elements from the following relationships: 


Na, Nb;Na;; — Nb,Na; 
Nb} = Nb; = Nb, is ies ts ble 


Nats Nats 
Nalz = Na = Na Nais - Na;jNats = Na,jN dis 
ij Uy cooN = N 
ats ats 


Ne} = Ne; — Ne, me 7 NcjNa,, — NcsNag; 
ats N ats 

5- We apply the stopping criterion of the direct Simplex algorithm on the objective function row. If the 
objective function is of the maximize type, the objective function row in the table must not contain any 
positive value. (But if the objective function is of the minimization type, the objective function row in the 
new table must not be contains any negative value), if the criterion is not met, we repeat the same steps 
until the stopping criterion is met and we obtain the desired ideal solution. 

Steps of the binary simplex algorithm: 


a. We write the two models in basic form by adding or subtracting additional variables or using synthetic 
variables and isolating the non-restricting variables. 
Basal form of the original model: 


Find 


NZ = NeyxXy + NegxX2 ++ + NCpXp + Ou, + Ouz +++: + Ou, — Max 
Constans: 
NQy1X1 + Nay2X2 + +++ + NaynX%n + Uy, = Nb, 
NQg4X1 + NQg2X2 +++ + NAgnX%n + Uz = Nb 
NQmiX1 + NAmn2X2 ++: + NAmyXn + Um = Nbm 
xj 20 ;j=1,2,..,n 
u; 20 ;i=1,2,...,m 
Here we do not require that Nb; = 0. 
Basic form of the dual model: 
Find 
NL = Nbyy, + Nboy2 +++ + NbiYm + OV, + OV2 +++ + Ov, — Min 
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Constans: 
Nay1y1 + Naz yz ++ + NOmi Ym — V1 = Ney 
Nay2¥1 + Nagzy2 +++ + NAm2Vm — V2 = Nez 


NaynY1 + NQgn¥2 + + + NOmnVm — Vn = Nn 
y, 20 ;i=1,2,..,m 
vy, 20 5j=1,2,...,n 
Here we do not require that Nc; = 0 
The two models have the same coefficients, and the matrix of instances of the dual model is the transpose of the 
matrix of instances of the original model. We write the two models in the following binary table: 


Table No. (3) Standard format for the original and companion models 


objective Max 
function Constan: 
constants col 


1 NaQy4X1 + Nay2X2 +++ NQynXn + Uy, = 
2 Nay4X, + Nay2X2 + mee + NQ27Xn + U2 


m NQm1X%1 + NAm2X2 +2 + NAmnXn + Um 


Non-negativity Hy) XQ, 22) Xp Uy, Ug, +) Um = 
constraints 


objective Min 
function Constan: 
constants col 


1 Nay1¥1 + NQg1¥2 +++ + NOmiVm — Vy a 


2 Nay2y1 + Nazz2y2 + +++ + NAm2Vm — V2 = 


n Naini + NQon¥2 ++ + NAmnY¥m — Yn = 
Non-negativity > 
constraints 
b. We place the variables and coefficients of the original model in the modified simplex table, and we place 
the variables of the dual model outside the table as follows: 


Table No. (4): The binary table for the original and dual models according to the modified Simplex algorithm 
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Binary simplex algorithm for the original and dual models: 


From the modified simplex algorithm of the original model, we obtain the optimal solution of the original model 
when all the elements are in the last row (the objective function row of the original model) Nc; < 0 ;j = 
1,2,...,n and at the same time all the elements are in the column The last (associated objective function column) 
Nb; = 0 ;i =1,2,...,m and we get the optimal solution for the dual model when all elements in the last column 
(associated objective function column) are Nb; => 0 ;i = 1,2,...,m and at the same time the last row (the objective 
function row of the original model) Nc; < 0 ;j = 1,2,...,n (because it will correspond to NG = —v;) which are 
the two conditions Same for both models. Therefore, when searching for the optimal solution for both models 
together, we must work to make all elements Nb; => 0 ;i = 1,2,...,m and to make all elements No S0;j= 
1,2,...,n to achieve this, we rely on one of the two models, put its variables and coefficients in a table, and place 
the dual model in an external frame of that table. In general, we find that the necessity of placing the two models 
in a short table does not allow us to get rid of the negative constants on the right side, and therefore the general 
case of the previous binary table can it must include negative constants Nb; < 0 , and the elements of the last row 
can include positive elements Nc; > 0 , so when searching for the optimal solution for the two models, we must 
work to address these elements based on one of the two models. 


Depending on the original model, we do this in two stages: 
The first stage: 


We make the constant Nb; non-negative, which corresponds to obtaining a non-negative basic solution for the 
original model. 


The second stage: 


We make all elements of the objective function row non-positive (in the case of the objective function, 
maximization), and this corresponds to obtaining the optimal solution required for the original model. 


Based on the dual model, we do this in two stages: 
The first stage: 


We must make the elements of the dual model objective function column Nb; > 0 ;i = 1,2,...,m ,The last row 
is non-negative 
The second stage: 


We must make the free constants for the dual model —Nc; non-positive, and this corresponds to obtaining the 
optimal solution for the dual model. We explain the above through the following example: 


Find the optimal solution for both the following neutrosophic linear model and its dual using the binary algorithm 
Example: 


Find: 
[5,8] x, + [3,6] x2 — Max 


Constans: 
2x, + 3x2 < [14,20] 


14 
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2x, +.x2 < [10,16] 


3x, < [12,18] 
3x, < [15,21] 


x, 20,x, 20 


We form the binary table of the model and the dual model: 


Table No. (5) The original model and its dual mode 


2X1 + 3X2 
2X1 + Xp 


3X, 


IA | IA 


IA 


IV} IA 


Non-negativity 
constraints 


constants 


1 


2y1 + 2y2 + 34 


3y1 + yz + 34 


Non-negativity 
constraints 


In the following table, we wrote the two models in standard form: 


2X, + 3x2 + Uy, 

2X1 + x2 + Uz 
3x2 + U3 
3x, + Uy, 

X%4,%X2, Uy, U2, Uz, Ug 


Non-negativity 
constraints 


constants 


1 


2y1 + 2y2 + 3Y4 — Vy 


3y1 + 2 + 3¥4 — V2 


15 
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Non-negativity > 
constraints 


Table No. (6) Standard format for the original model and the dual model 


We notice from the table that the standard form of the original model includes a ready-made base of additional 
variables u,,U2,U3,U4, but for the dual model there is no ready-made base. Therefore, we multiply the two 
restrictions by (-1) and we obtain the basic form of the dual model. 


The following table shows the basic form of the original and dual models: 


a 


2X, + 3x2 + uy, 


Table No. (7): The basic shape of the original model and the dual model 
2X, + X2 + U2 
3X2 + U3 


4a OCOCOCSOCSCi ay SCSCS~* 


Non-negativity X4,%X2, Uy, Uz, Uz, Uy 
constraints 
objective function 
enstants 
1 


Non-negativity 


= 


—2y, — 2y2 — 3y4 + Vy 


—3Y1 — V2 — 3¥4 + V2 


constraints 


We put the two models in the modified Simplex algorithm table and we get the following table: 
Table No. (8): The binary table for the original and dual models according to the modified Simplex algorithm 


Non-basic 
ibrable 
basic vibrable 
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Non-basic 
vibrable 
basic vibrable 


The first stage: 
1- For the original model: 
Since the values in the constant’s column are all positive, we study the values in the objective function row and 
determine the largest positive value. We find: 
max(([5,8], [3,6]) = [5,8] 


It is an expression of the variable x,. This means that it will enter the base. To determine the element that will exit 
from the base, we calculate the index N@, where: 


[14,20] [10,16] [15,21]] _ [15,21] 


N@ Emin a 3 


= [5,7] 


We find that the pivot column is the column of the non-base variable x,, meaning that the variable x, will enter 
the base instead of the variable u,, and the pivot element is the element resulting from the intersection of the pivot 
row and the pivot column, which is (3) 


We perform the switching between variables using a modified simplex algorithm. 


2- For the dual model: 


We study the elements of the objective function row. We notice that all the values are positive. Therefore, we 
study the elements of the constant’s column. We find that they are all negative values. We choose the most negative 
of them, which is (—[5,8]) which is the row of the base variable v,, so its row is the pivot row. To determine the 
pivot column and the pivot element, we calculate the index N@'where: 


[14,20] [10,16] [15,21]] _ [15,21] 


N@' € Max a | 3 


So, the column of the non-base variable u, is the pivot column, meaning that the variable u, will enter the base 
instead of the variable v,, and the pivot element is the element resulting from the intersection of the pivot row and 
the pivot column, which is (-3). We perform the switching between the variables using the modified simplex 
algorithm, from (1) and (2) We get the following double table: 


Table No. (9): The binary table for the first stage, the solution according to the original and dual models 


Non-basic vibrable 
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Non-basic vibrable 


The second phase: 
We apply the stopping criterion of the algorithm 
For the original model: 


Since the values in the constant’s column are all positive, we study the values in the objective function row. We 
notice that there is a positive value, which is [3,6], meaning that we have not yet reached the optimal solution. 
Therefore, we specify the pivot column, which is the column of the variable x2 corresponding to the only positive 
value in the objective function row. [3,6] In order to determine the pivot row and the pivot element, we calculate 
the index N@, where: 

[4,6] [0,4] [12,18]] [4,6] 


N@ € min |——,—_, 
min 3 1 


3 3 


It corresponds to the base element uw, so its row is the pivot row and the pivot element is (3). We swap between 
the variables using the modified simplex algorithm. 


For the dual model: 


We study the elements of the objective function row. We notice that all the values are positive. Therefore, we 
study the elements of the constant’s column. We find that there is a single negative value, which is (—[3,6]), which 
is the line of the base variable v3, so its row is the pivot row. To determine the pivot column and the pivot element, 
we calculate the index NO’ where: 

[4,6] [0,4] [12,18]] _ [4,6] 

a? = = 

So, the column of the non-base variable y, is the pivot column, meaning that the variable yy, will enter the base 
instead of the variable v2, and the pivot element is the element resulting from the intersection of the pivot row and 
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the pivot column, which is (—3). We perform the switching between the variables using the modified simplex 
algorithm, from (1) and (2) We get the following double table: 


Table No. (10): The binary algorithm table for the second stage 


Non-basic vibrable 


basic vibrable 


Non-basic vibrable 


We apply the stopping criterion of the algorithm: 


1- For the original model, we study the elements of the objective function row until the criterion for stopping 
the algorithm is met, which is the absence of any positive element 
2- For the dual model, we also study the elements of the constants column until the criterion for stopping 
the algorithm is met, which is the absence of any negative element 
3- We find that the criterion has been met and thus we have reached the optimal solution 
The optimal solution of the original model is: 


4 4 
xZE [5-2] ,U5 € [5-2] ,uz € [8,12] , xf € [5,7] ,uj =ug =0 
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The value of the objective function corresponds to: 


Z* = max Z € [29,68] 


The optimal solution of the dual model is: 
yi € [1.2] ,y4 € [1,6] v2 =y2 = y3 = 17 = 0 
The value of the objective function corresponds to: 
L* = MinL € [29,166] 
We note that: 
Z* = Max Z € [29,68] < L* = MinL € [29,166] 
This solution is acceptable according to the following theory: 


If (X4,X2,...,X,)is an acceptable solution to the original model of type Max _ and 
(V1), «++» ¥m)was an acceptable solution for the dual model of type Min, so the value of the objective function of 
the original model does not exceed the value of the objective function of the dual model for these two solutions, 
that is, it is 


n m 
j=1 i=1 


This is for all acceptable solutions for both models (including the optimal solution) 
5. Conclusion and Results: 


From the previous study, we arrived at a solution for the original and dual models at the same time, which are 
neutrosophic values from which we know the minimum and maximum profit that we can obtain, because the 
interpretation of the optimal solution for the original model is that it gives us the best production plan that makes 
the value of that production as large as possible, within Available capabilities. As for the optimal solution for the 
dual model, it gives us the best values for the prices of raw materials, which, if used without waste, will also give 
us the best production plan, and the result is the maximum profit. 
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